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We study renormalizable extensions of the standard model that contain vector-like 
fermions in a (single) complex representation of the standard model gauge group. 
There are 11 models where the vector-like fermions Yukawa couple to the standard 
model fermions via the Higgs field. These models do not introduce additional fine- 
tunings. They can lead to, and are constrained by, a number of different flavor- 
changing processes involving leptons and quarks, as well as direct searches. An 
interesting feature of the models with strongly interacting vector-like fermions is that 
constraints from neutral meson mixings (apart from CP violation in K° — K° mixing) 
are not sensitive to higher scales than other flavor-changing neutral-current processes. 
We identify order l/(47rAf) 2 (where M is the vector-like fermion mass) one-loop 
contributions to the coefficients of the four-quark operators for meson mixing, that 
are not suppressed by standard model quark masses and/or mixing angles. 


I. INTRODUCTION 

The driving force behind many of the extensions of the standard model (SM) has been 
the hierarchy puzzle. However, there may be reasons that the fine tuning of quadratically 
large contributions to the Higgs mass from very high momentum scales is acceptable (e.g., 
environmental reasons). This motivates the consideration of simple extensions of the stan¬ 
dard model that do not address the hierarchy puzzle and do not introduce any additional 
fine tunings of parameters. 

The SM provides no understanding why there are three generations of chiral fermions 
with masses that require weak symmetry breaking or why they only come in such simple 
representations of the gauge group. Vector-like fermions can be much heavier than the SM 
fermions, since their masses do not require weak symmetry breaking. For some choices of 
quantum numbers, such vector-like fermions can Yukawa couple to SM fermions. Given our 
lack of understanding of the number of generations and their quantum numbers, it seems 
worthwhile exploring the possibility that vector-like fermions exist. Of course, their masses 
may all be much larger than the weak scale, but it is not unreasonable that one such vector¬ 
like representation has a mass light enough that it can influence the next generation of 
flavor physics experiments. This motivates a study of extensions of the standard model with 
additional fermions that reside in a single vector-like representation of the SM gauge group 
with a mass M. (For some earlier studies of such models, see, e.g., Refs. [TUTU].) We focus 
on vector-like fermions that can have Yukawa couplings involving the Higgs field and the SM 
fermions. This allows them to influence flavor physics. There are several such models. These 
models have been studied in the literature previously, but here we consider all such models 
and compare the patterns of deviations from the SM they would give rise to. These models 
predict a very specific pattern for the contributions of beyond the standard model (BSM) 
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physics to A F = 1 flavor-changing neutral-current (FCNC) processes compared to A F = 2 
neutral meson mixings. For large masses M the dominant order 1/(47tM) 2 contribution to 
the coefficients of the four-quark operators responsible for neutral meson mixing arise at one 
loop, and are not suppressed by SM quark masses and/or weak mixing angles. We compute 
this contribution to meson mixing in the 7 models that involve strongly interacting vector-like 
fermions. In addition to considering A F = 1 and A F = 2 flavor-changing neutral-current 
processes, we also briefly discuss violations of lepton universality and unitarity of the CKM 
matrix that arises from the corrections to the W boson couplings in these models. 

The light neutrino masses play no role in our analysis so we take the SM to contain 
massless left-handed neutrinos. We are most interested in experiments that have reach in 
vector-like fermion mass of more than ~ 10 TeV when the Yukawa coupling constants of the 
vector-like fermions to the standard model fermions are around unity. For such masses, 
effects originating from the non-unitarity of the 3x3 CKM matrix can be neglected. 

If the new fermion is in a real representation it can have a Majorana mass term which 
does not require a helicity partner in the spectrum. In that case, unless the Majorana mass 
is very large or the Yukawa coupling of the new neutral fermion to the SM neutrinos is very 
small, a SM neutrino gets an unacceptably large mass. Hence we do not consider models 
where the additional fermions come in a real representation of the gauge group j/] 

ATLAS and CMS have searched for strongly interacting vector-like fermions contained in 
several of the models discussed below HM5]. The bounds from the Run 1 data constrain 
the masses of vector-like quarks typically at about the M > 800 GeV level. These constraints 
come from the pair production of the vector-like fermions via their gauge couplings, and are 
(essentially) independent of the Yukawa couplings, which are the focus of this paper. 

In Section II we describe the models, Section III explores the constraints from measure¬ 
ments both in the lepton and quark sectors, and Section IV contains our conclusions. 


II. THE MODELS 

There are 11 renormalizable models with vector-like fermions in complex representations 
of the standard model gauge group, where the vector-like fermions have renormalizable 
Yukawa couplings to the SM fermions through the Higgs doublet. 

Models I and II have SU(2) L singlet, E, and SU(2) L triplet, T e , vector-like fermions. 
Their SU(3)c x SU(2)l x U(1)y quantum numbers are (1,1,—1) and (1,3,—1). These 
vector-like fermions Yukawa couple through the Higgs field to the SM left-handed lepton 
doublets L\ 

Model III has vector-like fermions, £, with same gauge quantum numbers as the SM 
left-handed doublets, (1,2,—1/2). They Yukawa couple to the SM right-handed charged 
lepton fields e R■ Model IV contains vector-like fermions, l, with (1, 2, —3/2) gauge quantum 
numbers, and they also Yukawa couple to e l R . 

Models V and VI contain vector-like fermions D and U with the same quantum numbers 
as the SM right-handed down and up-type quarks, (3,1, —1/3) and (3,1, 2/3). They Yukawa 
couple to the SM left-handed quark doublets Q l L . Models VII and VIII are similar with 
vector-like SU(2)l triplet fermions Td and T u with gauge quantum numbers (3, 3, —1/3) and 
(3,1, 2/3) respectively. 


1 There are two cases of this type. New fermions with SM quantum numbers (1,1, 0) and (1,3, 0). 
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There are three more models where the vector-like fermions Yukawa couple to the SM 
right-handed up-type, u l R and down-type d l R quarks. In Model IX the vector-like fermions 
q u d have the same quantum numbers as the SM left-handed quark doublets, (3, 2,1/6). In 
this case the new vector-like quark doublet couples to both u l n and d l R . Model X has a 
vector-like quark doublet q u with quantum numbers (3, 2, 7/6) that Yukawa couple through 
the Higgs doublet to the right-handed up-type quarks. Model XI has a vector-like quark 
doublet qd with quantum numbers (3,2,—5/6) that Yukawa couples to the right-handed 
down-type quarks. 


A. Lagrangians 


The new terms in the Lagrange density which are added to the SM Lagrange density are: 


Model I (1,1, -1) : £ B sm = E (tip - M) E - (\E R rfU L + h.c) , (2.1) 

Model II (1, 3, -1) : £ BS m = Tr [f e (tip - M) T e ] - (\H^T eR U L + h.c.) , (2.2) 

Model III (1,2, -1/2) : £ BSM = M) 1 - (A £ r HH l + h.c) , (2.3) 

Model IV (1, 2, -3/2) : £ BSM = Z (0 - Af) Z - (X,e l R H T e l L + h.c.) , (2.4) 

Model V (3,1, -1/3) : £ BSM = D (tip - M) D - (X 2 D R H^Q l L + h.c.) , (2.5) 

Model VI (3,1, 2/3) : £ BSM = U (tip - M) U - (A iU R H T e Q\ + h.c.) , (2.6) 

Model VII (3, 3, -1/3) : £ bsm = Tr [f d (tip - M) T d ] - + h.c.) , (2.7) 

Model VIII (3, 3, 2/3) : £ BSM = Tr [f u (tip - M) T u ] - (X 2 H T e T uR Q\ + h.c) , (2.8) 

Model IX (3, 2,1/6) : £ BS m = q u d (0 - M) q ud - (x\ u) u R H T e q udL + Xf ] d l R Hp[ udL + h.c), 

(2.9) 

Model X (3, 2, 7/6) : £ BSM = q u (0 - M) q u - (X t u l R H' f q uL + h.c.) , (2.10) 

Model XI (3, 2, -5/6) : £ BS m = Id (0 - M) q d - (Xid l R H T eq dL + h.c) . (2.11) 


Here the 577(2) triplets T e , T u and T d are represented by the two-by-two matrices 


( T- i/V 2 T 0 A 

V T_ 2 -T^/V2)' 

T- 1 / 3 /v/2 T 2 /3 A 

T-4/3 —T- 1 / 3 /V2 ) 

T = ( ^" 2 / 3 /02 T 5/3 A 

“ V T_ 1/3 -T 2/3 /y/2j' 


(T-xN 2 T_2 A 

V T 0 -T.i/v^;’ 

T-1/3/ a/2 _T_ 4/3 A 

T 2/3 —-T- 1 / 3 // ’ 

r = (T2/3/V2 T- 1/ 3 A 

u V t 5/3 -t 2/3 /V2J' 


( 2 . 12 ) 

(2.13) 

(2.14) 


where the subscripts denote the charges of the fermions. 

It is important to note that there are no new one-loop contributions to the ordinary 
lepton and quark mass matrices of the form 


"Vo- ~ 


A*A,u 

167T 2 


(2.15) 


where v ~ 174 GeV is the Higgs vacuum expectation value. This is due to the approximate 
chiral symmetry, e.g., e l R —» e lOL e l R in Model I, which is only broken by the SM Yukawa 
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couplings. This prevents such contributions to the light fermion mass matrix. So the A* 
couplings can be order unity without the need for a cancellation in the SM fermion mass 
matrix. 

The BSM effects caused by the vector-like fermions vanish as their mass M —>• oo. We 
are interested in very large M, greater than 10 TeV, and it is flavor-changing neutral- 
current processes that are sensitive to such BSM physics. In the SM, flavor-changing neutral 
currents are suppressed by small weak mixing angles and small SM quark or lepton masses. 
Furthermore, they do not occur at tree level. Since the neutrino masses are very small, 
charged lepton flavor violation in the SM is negligibly small. In the quark sector, despite their 
suppression, flavor-changing neutral-current processes that change flavor by one (A F = 1, 
e.g., K + —y n+vv) and by two (A F = 2, e.g., K° — K° mixing) have been observed. 

In the limit where the SM Yukawa couplings vanish, the SM has a U(3)q x U(3) Ur x 
U(3)d a X U{3)l x U{3) £r flavor symmetry. This symmetry forbids flavor-changing neutral 
currents. However in Models I-XI, even when the SM quark and lepton Yukawas are zero, the 
terms in the Lagrange density proportional to the Yukawa couplings Aj break that symmetry. 
That pattern of flavor symmetry breaking can be characterized using the spurion method. 
For example in Model V the BSM terms in the Lagrange density involving X t are not invariant 
under the U(3)q transformations Q l L —> V{Qi)ij Ql (repeated flavor indices are summed). 
However, if the A* also transform as A* —>■ H(Q/-J* ? A J , then the U(3)q symmetry is restored. 
This means that in Model V one-loop BSM physics associated with the high mass scale M 
generates meson mixing through a term in the effective Lagrangian of the form, 

•c£L„ ~ | (Qh,Qi) (Q'LYQ'l) , (2.i6) 


and similarly for the other models. These corrections are not suppressed by small quark 
masses and/or mixing angles, and we have not found expressions for them in the literature. 
(We take the BSM Yukawa couplings A; to be of order unity.) They are computed for 
Models V-XI in Sec. IIIB4 Corrections to meson mixing of order \\*/{AnM) 2 that are 
also suppressed by weak mixing angles and/or quark masses were considered for Model V 
in Ref. H6l. 


B. New interactions with gauge bosons 

The new BSM contribution to the Z coupling and to the W arise at tree level, after 
integrating out the heavy fermion(s). They can be obtained ether by calculating Feynman 
diagrams or by diagonalizing the 4x4 fermion mass matrices as was discussed in Ref. m 
Here we explicitly show how to obtain the BSM contributions to the Z couplings for Model I. 
(The same method can be applied to the other models.) 

After electroweak symmetry breaking, the mass matrix in charged lepton sector becomes, 

An® = -eiM e AB en + h.c. , (2.17) 

where roman capital indices A and B go over {0,1, 2, 3}, e is defined by e = (E , e 1 , e 2 , e 3 ) T , 
and 

/MO 0 0 \ 

A^u m e i 0 0 

X^v 0 m e 2 0 

^ \* 3 v 0 0 m e3 ) 


M 


(2.18) 
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Here we assumed, without loss of generality, that the charged lepton fields e l L (i.e., the lower 
components of the doublets L 1 l ) are eigenstates of the charged lepton mass matrix in the SM 
(i.e., in the A* —> 0 limit). This matrix is diagonalized by the 4x4 unitary transformations 

v L Ri 

&'l,r = Vl,r£l,r , (2-19) 

where the prime denotes a mass eigenstate field. So 


V£M 6 Vj? = M 


l diag 


Up to corrections suppressed by (■ v/M) z 


V,f 


( 1 

-A \v/M 
—\* 2 v/M 
\ -A %v/M 


Ai v/M X 2 v/M A 3 v/M^ 

1 0 0 

0 1 0 

0 0 1 ) 



*-4x4 • 


( 2 . 20 ) 


( 2 . 21 ) 


Consequently the masses of the charged leptons in the SM part of the Lagrangian are 
approximately equal to the charged lepton masses, and the vector mass parameter M is 
approximately the heavy vector-like lepton mass. 

The Z boson coupling is, 




—gzZ^ 



(Vi)A0 7“K 


)oB z'l 



( 2 . 22 ) 


where the ellipses denote terms not containing the matrix V£. Here gz = \J g\ + g 2 and gi ,2 
are the gauge couplings of U(l)y and SU(2) L , respectively. There is some ambiguity in how 
the terms are organized since V[ is unitary. We have written the Z couplings involving V£ 
so that the order (v/M) 2 corrections to the SM Z couplings to the charged leptons can be 
read off using the order (v/M) terms in V£ that we have explicitly calculated. From this it 
follows (removing the hats and primes) that the BSM couplings of the Z boson to the light 
mass eigenstate charged leptons are, 

Model I : £<f s > M - - W , (2.23) 


with rriz being the Z boson mass. For simplicity, hereafter we remove the primes used to 
specify the mass eigenstate fields. 

In the same way, the new Z couplings to the quarks and charged leptons are obtained for 
the other models. (The 4x4 unitary matrices are given in Appendix |Aj) The results are 


Model II : 

Model III : 

Model IV : 



(2.24) 

(2.25) 

(2.26) 
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Note that in Models III and IV the BSM tree-level Z couplings differ only by an overall sign. 
Similarly for hadronic models, 


Model V : £g> M = - X ( , 


hJ 


Model VI : £g> M = X 




5# 2 


A * V4 , , , v z 

a Ll a L J 


fcM 2 


Model IX : £g> M = - X 


A-^A• X\ a> \) a> *mi _. 

J u' R ru> R - * ' z <T I tf‘<P R \Z (2.29) 


(d) \ (d)* 2 




Model X : £g> M = X 


h3 


A,A*m| 

w 


gzM 2 Hl H g z M 2 
^ 1 rY u r 1 


Model XI : £g> M = - £ 


A*A*m| 


^ • 


(2.27) 

(2.28) 


(2.30) 

(2.31) 


So far, we have worked in the mass eigenstate basis most convenient for expressing the 
Z couplings. For example, in Model V we have worked in a basis where in the Aj —>■ 0 limit 
the down-type quarks in the Lagrangian are mass eigenstates, while in Model VI we have 
worked in the basis where the down-type quarks in the Lagrangian are mass eigenstate. In 
Model IX the couplings A ^ correspond to the down quark mass eigenstate basis and the 
couplings A < -“' ) are in the up-quark mass eigenstate basis. 

In case of vector-like triplets, i.e., Models VII and VIII, it is not possible to choose a basis 
where the CKM matrix, V, is absent from the Z couplings]^] Choosing the up-type quarks 
in the Lagrange density to be mass eigenstates (in the Aj —> 0 limit) implies for Model VII, 


Model VII : £g> M - - V 


g z M 2 


sj/rX + b E 4 vl r V jn dl)z „ . (2.32) 


Similarly, in Model VIII the result is 
<z) _ 


Model VIII : £g' M = V 


^ 0 -) (\sird + X r<r L )z„. (2.33) 


We can transform these Z couplings to the basis where the down-type quarks are mass 
eigenstates in the A* —» 0 limit by redefining the couplings A, —> A kV ki . In this basis 

Eqs. (2.32) and (2.33) become 


Model VII : £g> M = - X ( 


hJ 


X «” V mi Y v] n ul + - cf L Yd{ z„ , (2.34) 


v m,n 


and 


Model VIH:£g> M = -X(^g? 


hO 


X 5? V mi Y V] n ul + S l YA\ Z „. (2.35) 


2 More precisely, the unitary matrix V is the CKM matrix in the Xi —> 0 limit. 
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There are corrections to the W boson couplings as well. We write the couplings of the W- 
bosons to the left-handed light lepton and quark mass eigenstates in terms of 3 x 3 matrices 
X and Y a^J 


£(W) 


92 

V 2 


K [PirXaei] 


■^W+[u‘ L rV,,d’ L ]+h.c.. 


(2.36) 


In the SM limit, A" is the inverse of the PMNS matrix and Y is the CKM matrix. Because of 
the mixing with the vector-like leptons in the extensions of the SM discussed in this paper, 
X and Y are no longer unitary matrices. It is straightforward to express X and Y in terms 
of the components of the 4x4 diagonalization matrices in the various models. 

In this paper we focus on violations of lepton universality and violations of unitarity of 
the CKM matrix. For that purpose we compute the quantities. 


** = (**’)»= E x * x h’ 

3 = 1 , 2,3 

&» = WL = E W,. 

3= 1,2,3 


(2.37) 

(2.38) 


at quadratic order in ( v 2 /M 2 ), neglecting terms of higher order. We find in the models with 
vector-like leptons that 


Model I : 

Model II : 
Models III, IV : 


Rk ^ 1 - |A fe 



~ 1 . 


M 2 ’ 

|Afc|" + |A' fc | 2 - \* k \' k - AfcA^ 



(2.39) 

(2.40) 

(2.41) 


In Model II we have worked in a basis where the charged leptons e J are mass eigenstates in 
the SM part of the charged lepton mass matrix and so A'- = JA X,Uij, where U is the PMNS 
matrix. 

In the models with vector-like quarks (neglecting the off diagonal elements of the CKM 
matrix) 


Models V, VI : S lm ~ 5 lm - , 

(2.42) 

Models VII, VIII : S lm ~ S lm + ^A ; *A m ^ , 

(2.43) 

Models IX, X, XI: S lm ~ 5 lm . 

(2.44) 


At this order in the v 2 /M 2 expansion, X) m = (Y^Y)i m = S) m . 


III. EXPERIMENTAL CONSTRAINTS 

For leptonic-extension models, i.e., Models I - IV, the constraints from fi —» e conversion, 
fi y 3e and r —» 3e are important. They are induced by the tree-level Z couplings, which 


3 There are also right-handed currents in some of the models. 



we derived in the previous section. On the other hand, for hadronic-extension models, i.e., 
Models V - XI, meson mixing, such as K°-K° and D°-D°, are induced at one-loop level. 
In addition, A F = 1 FCNC processes, such as K —* 7 tuu, K —* /i + /i _ , B s —> i + i ~, etc., 
are induced at tree level. We will derive effective Hamiltonians which are relevant for these 
processes and discuss the current experimental bounds, as well as future prospects. (For 
some recent studies of constraints on some of these models, see, e.g., Refs. pT8ld2C7j .) 

For upper bounds on lepton flavor violating processes we quote bounds at 90% CL, as 
do most experiments, whereas for other measurements we quote the la limits that are 
approximately 84% CL as one-sided bounds for Gaussian distributions. 


A. Leptonic models 


The flavor violating tree-level couplings in Eqs. (2.23)—(2.26) give BSM tree-level con¬ 


tributions to the flavor-changing neural-current processes; for example, in the muon sector, 
to the /r —)■ 3e rate and muon conversion to an electron in the presence of a nucleus. The 
amplitudes for these processes are both proportional to the same combination of parameters 
A 1 A 2 /M 2 . The flavor diagonal terms in Eqs. (2.23)—(2.26) give rise to violations of univer¬ 


sality in lepton couplings to the Z. For typical couplings A these are not as sensitive to large 
values of M as the flavor-changing charged lepton neutral currents are. Note that in these 
models the radiative decay /d —> ey does not arise at tree level, only at one loop. 

Let us first discuss // —» e conversion in the presence of a nucleus. It is among the most 
powerful probes of charged lepton flavor violation beyond the standard model. The 
conversion rate to an e~ in the presence of a nucleus N is usually quoted as a branching 
ratio normalized to the SM weak interaction, /d~N —>• u^N' capture rate. At present the 
most stringent bound is Br(/i — > e conv. in Au) < 7 x 10 -13 at the 90% C.L. [21]. Future 
experiments will have a dramatically improved sensitivity to a branching ratio Br(/d —> 
e conv. in Al) ~ 10" 16 [221 123] • 

The reach, in mass scale for new physics, of the next generation charged lepton flavor 
violation (CLFV) experiments that search for /i to e conversion on Al and for the radiative 
decay jd —» ey was discussed in Ref. [24]. They assume an effective Lagrangian of the form, 


Cr 


- tpllv 

CLFV = , 1 N A2 eL + 


K 


and present the reach of these experiments in the k — A plane. For k 1 the limit on A from 
conversion on gold is close to 1 x 10 3 TeV while future planned experiments for conversion 
on aluminum are sensitive to A ~ 7.2 x 10 3 TeV. In Models I - IV the weak radiative 
decay proceeds at the one-loop level and so indeed k 1. Taking into account left-handed 
and right-handed quark currents which couple to Z boson, the relations between A and the 
vector-like lepton mass M in Models a = I, II, III, and IV are 


M (a) 


A 


44 , 


IA2 Ai 

3 


(2 sin 2 9 W — 




1 

2 



(3.2) 


with r)f,he = = 'rjjiXe = 1 and Z, A are atomic number, mass number of a 

nucleus, respectively. To obtain this expression we used the vector part in the quark current, 
neglecting the axial current. For Models III and IV, the right-handed lepton current occurs 
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but it gives the same contribution to the rate as the left-handed one. On gold and aluminum 


Eq. (3.2) implies that, 


A/fi = 0.31A ^“-LlAaAil, M$ = 0.28A yjrjj*Xe |A^|, (3.3) 

We obtain from Fig. 2 of Ref. m that a future branching ratio limit —> e conv. in Al) < 

10 -16 will imply A > 7.2 x 10 3 TeV for large n, yielding M j \Jr^Xe |A 2 Ai| > 2.0 x 10 3 TeV. 
The current limit on p —>■ e conversion in Au implies that A > 9.5 x 10 2 TeV (24], which 


according to Eq. (3.3) implies the limits, 

Af<“> > x 2.9 x 10 2 TeV. 


(3.4) 


Charged lepton flavor violation is also constrained by the muon decay p —> 3e and 
similar r decays. In the models a = {I, II, III, IV} the rates for these processes normalized 
to a leptonic weak decay mode are, 


and 


T(e* -+ e J e J e J ) 
T(e* —> v'eve) 

T(e* -+ e j e k e k ) 


(a) 

m 

sg 2 f 


I Aj Xj | 
M 4 


(a) 

8 G 2 f 


| AjAj 


T(e* -> z+eP e ) 8GJ, M 4 ' 

For the decays where the leptons in the final state are all of the same type, 

„(ni) 


V? = A Vi 1] = 24 + Kr , 


Vi 


= nl + 2k 2 


R i 


(3.5) 


(3.6) 


(3.7) 


where kl = —1/2 + sin 2 6w and kr = sin 2 9\y, and 9w is the Weinberg angle. Numerically: 
rj[ ]> ~ 0.20, ~ 0.049 and r/[ T11 ^ = i][ lx ^ ~ 0.18. For the decays where two types of leptons 

occur in the final state, 


(I) (III) (IV) 2 , 2 

V2 = V2 =V2 =K L + K R , 


(II) 1 (I) 
V2 = 4*72 


(3.8) 


Numerically; 77 ® = r /.® 1 ' 1 = r/®"* ~ 0.13 and 77 ®^ ~ 0.031. The 90% CL experimental limit, 
Br(/i —>■ 3e) < 1.0 x 10 ~ 12 [23| implies in Model I that M/-J\ AiA 2 1 > 1.2 x 10 2 TeV. While 
the limits Br(r —> 3/i) < 2.1 x 10 -8 and Br(r —)• 3e) < 2.7 x 10 -8 [26] imply in Model I that 
M/ a/|A 2 A 3 | > 6.2TeV and M/ y 7 1A 1 A 3 > 5.8TeV. The limits on the rates for the r decay 
channels r —» p,ee and r —> e/i /2 give slightly weaker limits because the r/ 2 ; s are smaller than 
the 77 ! ’s. 

The Mu3e experiment expects to reach in the absence of a signal the 90% CL limit, 
Br(/i —> 3e) < 4 x 10 -16 (27], yielding in Model I, M/ y/|AiA 2 | > 8.2 x 10 2 TeV. The Belle II 
sensitivities for r decays to three charged leptons are estimated at the few times 10~ 10 level, 
and several channels will give comparable sensitivity. The bounds Br(r —y 3e) < 4 x 10 -10 
and Br(r —> 3 fi) < 4 x 10 ~ 10 [2H] would yield in Model I, M/yJ |Ai j2 A 3 | > 17TeV. 

Competitive constraints to the above processes also arise from the upper bounds on the 
r —> en, r —>■ /nr, r —> ep, and r —>• pp branching ratios. The vector-like fermions generate 


T(r -+ e*7r) 


(a) |A 3 A / 2 ml 2 
/3 M A 256 7r ^ ’ 


(3.9) 
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where % — 1 , 2 , 77 ® = 4773 ^ = ^g 1111 = 77 ^ ) = 1 , that is, 773 ^ = t)^/(k 2 + /c^), and we 
neglected 77 ?. 2 /m 2 . Similarly, 


T(t —y e l p) 


(a) |A 3 Ai| 2 (ml 
% M* - 


m 2 p ) 2 (ml + 2m 2 p ) 
256 n ml 


(l - 2 sin 2 6 >w) 



(3.10) 


Of the current 90% CL experimental limits [2Ri I3T] the strongest bounds arise from Br(r —y 
en) < 2.2 x 10~ 8 , which implies M/<J |AiA 3 | > 7.0 TeV in Model I, and Br(r —» pp) < 
1.2 x 10 -8 , which implies M/ a/|A 2 A 3 | > 7.4TeV in Model I. The bounds from r —* ep and 
r —y pn are only slightly weaker. 

Concerning future sensitivity, the expected Belle II limits are [28]. 

Br(r ->• evr) < 4.0 x 10 ~ 10 , Br(r ->• pn) < 4.9 x 10 ~ 10 , 

Br(r —>■ ep) < 3.1 x 10 ~ 10 , Br(r ->• pp) < 2.0 x 10 ’ 10 . (3.11) 


These imply that the strongest expected bounds in Model I will be M/ a^/| Ai A 3 1 > 19 TeV 
from r —» en, and M/sJ |A 2 A 3 1 > 21 TeV from r —>• pp. The expected sensitivities of the 
other channels are only slightly weaker, thus we can have high confidence in the experimental 
reach, but not in which channel will give the best bounds. 

For generic BSM Yukawa couplings, A, the reach for new physics is much greater for exper¬ 
iments that search for charged lepton flavor violation than those that seek violation of lepton 
universality. However, it is possible for non-generic A’s that the violations of universality are 
more important. We close this section by briefly commenting on the implications of the pow¬ 
erful constraint on e-p charged current universality coming from pion decay 7 r + — y e + u('y), 
p + v(^). The latest experimental result gives g e /g p = 0.9996 ± 0.0012 [32j where are the 
charged-current couplings of e and /id] In Models I - IV, g e /g^ = R\/R- 2 , which implies, 
for example, in Model I that, 


M > y/| |A 2 | 2 - | Ax I 2 1 x 4.4 TeV . (3.12) 

It should be possible to improve this bound by a factor of two in the future. Tau decay, 
on the other hand, gives constraints on other charged-current couplings, such as g T /g^ = 
1.0011 ± 0.0015 |3Tj, which yields, e.g., in Model I, 

M > ^/||A 2 | 2 - |A 3 | 2 | X 6.2 TeV. (3.13) 


B. Hadronic models 

To constrain Models V - XI, FCNC processes in the quark sector are most important. 
(For a recent study of Z-mediated FCNC effects, see, e.g., Ref. [35].) We focus on lep- 
tonic and semilcptonic decays and neutral meson mixing. In these models, constraints from 
nonleptonic decays are weaker. 


4 The constraint from Kaon decay JM] is weaker. 
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1. Meson decays involving a vv pair 


For kaon decays involving a neutrino-antineutrino pair, such as K + —y i r + zzzq the effective 
Hamiltonian is 


V = c„ ? (s L 'y ll d L ){i'i L Yi'i L ) + <4 ) j 


(3.14) 


Cy,V 

~Y 


^2(s't„d)(v iL Yv i L) + ~y , 


(3.15) 


where in the second line of Eq. ( |3. 14 ) the subscripts V, A on the coefficients refer to the 
fact that these coefficients are for the vector and axial quark currents and c v y = d v + c„, 
c u ,a = d v — c v . (See Appendix [B]) Only the vector part of the quark current contributes to 
this process. In the SM d v = 0 and at NNLO, |c£ S ^| = |c]> SM! | ~ 1/(89TeV) 2 , with a few 
percent uncertainty [3&Hl0j. The hadronic models modify c„y additively. Obviously in this 
case it is convenient to work in the basis where the down-type quarks are mass eigenstates]^] 
The tree-level Z exchange BSM contributions to the coefficient c u y in models a = {V, VII, 
VIII, IX, XI} have magnitude 


-v,V 


(a) IVA2I 

- ’fej yy 


(3.16) 


where for the FCNC down-type BSM contributions 


47 = 24V = 4T = 4? = 4? = 1. (3.17) 

Note that for Model IX it is the A^’s that occur in this expression. Models VI and X do 
not contribute to this process through tree-level Z exchange. 

The uncertainty of the SM prediction, Br (SM ^(/l + —» n + ui i ) = (7.8 ± 0.8) x lO^ 11 , is 
dominated by that of c|f^, and not by its relation to the measured rate. Hence, the 

measurement Br(/t + —>■ n + uu) = (1.7±1.1) xlO -10 [HJ implies 0.8 < |l+c^ a y/c[f^| 2 < 3.6, 6 * 
In Fi g, [l] the bl ue-shaded region shows the la allowed region in the parameter space or 

M/\Jriy, | A 1 A 2 1 and 8 = arg {c^y/. If the BSM contribution is aligned (constructive 
interference) with the SM, we find 


M M > yi?zJ |AiA 2 x 66TeV. 


(3.18) 


This is not a bound in the strict sense. Anti-alignment, 8 ~ 7r, is possible for lower values of 
M. But it is more indicative of the new physics reach, since a larger region in 8 is allowed 
for values of M satisfying Eq. (3.18). 

The upcoming 10% measurement of Bt(K + —> 71 + vv), if consistent with the SM, would 
result in the (la) constraint 0.9 < |l + c'"y /c[ S y 1 ' 1 1 2 < 1.1, improving the sensitivity to BSM 
physics to about 2.8 x 10 2 TeV. This is shown as the red-shaded region in Fig. [l] Again 
we emphasize, as is clear from Fig. [lj that quoting the weakest possible bound would not 


5 Going forward, for different processes we work in whatever basis makes the theoretical expressions simplest. 
This means that the A parameters are not always the same, but are linearly related through the CKM 
matrix. 

6 Despite the stated uncertainty, the probability that all 7 observed events were due to background was 

quoted as 0.001 0U- 
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FIG. 1. Constraints from K + —> on vector-like fermions in the M/\J |AiA 2 1 vs. <5 = 

arg (c^y/c^ S ^). The currently allowed la region is blue-shaded, corresponding to 0.8 < |l + 
cH/c^i 2 < 3.6, whereas a future 10% measurement in agreement with the SM, 0.9 < |l + 
c lt/ c fr| 2 < 1-1) would constrain M and 5 to the red-shaded region. 



give the most useful impression of the scale sensitivity, as a solution around cj/y = —2c^p 
always remains, corresponding to a relatively low new physics scale. A 3 cr deviation from 

the SM is possible for M*®) j ^rj^l | AiA 2 1 < 1.7 x 10 2 TeV. 

For the similar 3rd-2nd generation transition mediated by b —y suv, the SM prediction 

) by d —> s and s —> b replacements 
is |c^l — 1/(9 .8 TeV ) 2 [02]. This process has not been observed yet, and the best current 
bound, Br(T> —> KuD) < 1.6 x 10 -5 , is about 4 times the SM prediction j03j. This yields 

for new physics aligned with the SM contribution, M ^ r/f/j |A 2 A 3 1 x 6.9 TeV. In case 
the new physics is anti-aligned, interfering destructively with the SM contribution, M^ > 

\J 77 ^] |A 2 A 3 1 x 4.0 TeV. If the rate is at the SM level, Belle II expects to measure Br(£> —» 
K*uu) with about 30% uncertain ty [28lj Q in creasing the probed mass scales to about 19 TeV 

(a 3a signal possible for M (a> /\Jr /S/j |A 2 A 3 | <11 TeV). 

For B —> nuu, Belle II expects to reach a sensitivity at the 1 x 10 -5 level [25], which will 
provide weaker bounds than B — » 7 t£ + £~ in the models considered in this paper. 


for the coefficient of the operator obtained from Eq. (3.14 


' At the SM level, B —> Kvv is expected to get large backgrounds from B —> K*vv [28]. While Br(I3 —> 
Kvv) depends only on C +v , Br(B —► K*vv) also depends on C- v , slightly complicating the analysis. 
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2. Meson decays to an t + l pair 


The effective Hamiltonian differs from Eq. (3.14) in that the coupling to the £ + £ pair 


can either be left- or right-handed. We write the effective Hamiltonian as 

TL = (cj Qi + c- Q [), 

i=9,10 

where up to normalization the conventional choice of operator basis is 

Q 9 = (b L ^s L ) ihn ), = (bia,s R ) (hn ), 

Qio = ( b L ‘J^s L )(£'f'y 5 £), Q' w = (b R 'y ll s R )(l'f'yst ), 


(3.19) 


(3.20) 


with obvious replacements for b —y d or s —y d decays. We use the notation q and c' 
to emphasize that these are dimensionful couplings, containing all terms multiplying the 


four-fermion operators Qi in Eq. (|3.20|). 
For B s —y n + fj,~ 

—y fi 


the SM gives (c^ 


~ l/(17TeV) 2 and M ' 1 1 ~ 0. The amplitude for 


B s —y is proportional to Cio,a = Cio — c' 10 . In the models with vector-like fermions 

the BSM contributions to the coefficients cio,a in models a = {V, VII, VIII, IX, XI} have 
magnitude 


J a ) 


(a) 


C 10,V — Vzd 


|A 2 A 31 

AM 2 


(3.21) 


and the r]zd coefficients are the same as in Eq. (3.17). The LHCb-CMS combination of their 
measurements, Br(H s —y /r + /r~) = (2. 8-ae) x 10 -9 , is quoted as the SM prediction times 
0.76^o)i8 [4I|- The la range of the measured B s —y /r + /i~ rate is slightly outside the SM and 

corresponds to the region, 0.96 > |1 + | 2 > 0.58. The allowed one and two sigma 

regions of BSM parameter space are plotted in Fig. fusing variables very similar to Fig. [lj 
If the BSM contribution is destructive (at la), 


M^> V / 4 a j|A 2 A 3 |xl8TeV. 


(3.22) 


A future 10% measurement 


would increase this sensitivity to ~ 40 TeV. Again we 
note that Eq. (3.22) is not a true bound on M but rather is meant to give a feeling for the 


reach in M of current experimental data on this decay mode. There is a tuned region of 
BSM parameter space near c^' ) _ 4 /c^ 1 y ) = —2 that corresponds to a lower value of M. 

For B d — y /i + /r _ the SM prediction is c,^ ~ 1/(37TeV) 2 . The LHCb-CMS combina¬ 
tion, Br (B d —y /r + /i~) = (3.91};®) x 10~ 10 is quoted as the SM prediction times 3.71};® 

If the new contribution is constructive to the SM, then 


M (a) > yjrizd | AiA 3 1 x 16 TeV . 


(3.23) 


Accommodating the current central value would require /yr/E/j |A 3 Ai| ~ 19 TeV. In 
the HL-LHC era, this measurement will reach an uncertainty around 20% of the SM predic¬ 
tion [45, 06], increasing the mass reach to above 60 TeV. 

The extraction of the short-distance part of the measured rate Br(A% —y /i + /r _ ) = (6.84± 
0.11) x 10~ 9 [47j is subject to considerable uncertainties. The estimate Br (K L —y /i + /r _ )sD < 
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FIG. 2. Constraints from B s — > pL + p>~ on vector-like fermions in the M /y r/^j | A2A3I vs. 5 parameter 
space. The angle delta is defined in a similar way as in Fig. [l| The blue-shaded region is allowed 
at lcr (indicating the mild tension with the SM), and the green-shaded region shows 2cr. 


2.5 x 10 9 |48j is about 3 times the SM short-distance rate which follows from |cf ^4 
1/(1.8 x 10 2 TeV) 2 . If the BSM component is aligned with the SM, we find 


M (a) > 



|Re(AiA 2 )| x 100 TeV. 


(3.24) 


Unlike the bound from K + —> 7 t + z/R the prospect of improving this is not good, and the 
uncertainties are greater. 

Recently LHCb established a strong bound Br(H° —> /i + /r _ ) < 6.2 x 10~ 9 at 90% CL [49], 
which is well above the SM level. In Models VI-X, 


Br (D° /i+p, ) = t d 0 flm D ml^l-Aml/m 2 D , (3.25) 


where 


u (VI) 

r iZu 


„(vn) = 2?? ( vm ) = J lx) = n (X) = 1 

tZu ^ IZu 'IZu 'IZu ^ ' 


(3.26) 


For Models VII and VIII we used the basis in Eqs. ( |2.32[ ) and (J2433J). Using f D = 
0.209 GeV [HU] , we obtain 


M (a) > \jrjf u lAiAal x 3.9TeV. (3.27) 

It is possible that in the HL-LHC era the experimental bound will improve by a factor of 

~ 20 [ 3 ]. 


3. Semileptonic decays to £ + £ pairs 

The bounds on K L —>■ n°£ + £~ are about an order of magnitude above the SM expectation, 
so the resulting constraints are weaker than those obtained from K —> uvu and K L —>• 














15 


LHCb recently measured Br(R + —* 7r + /i + /i _ ) = (2.3 ± 0.6) x 10 -8 [52], consistent with 
the SM prediction quoted as (2.0 ± 0.2) x 10 -8 . This was the first FCNC b —> d decay 
observed (other than mixing). Requiring that the SM rate is not enhanced by more than 
50% by new physics, using Appendix [Bj we find 

M (a) > |AiA 3 | x 30 TeV. (3.28) 


This bound is slightly stronger that that from Bd —> /i + /i _ , moreover, the current B + —> 
7 r + /r + /i _ measurement only used 0.9/fb data. As the measurement gets more precise, a ded¬ 
icated analysis of B + —> 7 r + /i + /i _ , possibly considering [dT(R + —> 7 r + /i + /x _ )/dg 2 ]/[dT(R —$■ 
7 iPv)/dq 2 } to reduce theoretical uncertainties, is warranted. 

The analysis of B —> K/i + fj,~ is very similar to B + —> 7r + /i + /i _ , while B —> K* /i + fi~ is 
more complicated and does not give better bounds. Averaged over l = e, /i, HFAG quotes 
Bt(B — > K£ + £~) = (4.8 ± 0.4) x 10 -7 pH]. The experimental uncertainty is smaller than 
the theoretical one (due to the form factors). Demanding a less than 30% modification of 

the SM rate, we find M^ |A 2 A 3 1 x 17 TeV. This happens to be very close to the 

bound in Eq. (3.22), but the prospects of improving that are better. 

The inclusive decay rate Br(R —>■ X s £ + £~) = (5.0±0.6) x 10 -6 [3T] depends in the models 
we consider on |C 9 1 2 ± |cg ( 2 and |C 10 1 2 ± |c/ 10 | 2 . In Models IX and XI, there is no interference 
between the SM and the new physics contributions, so the constraints are weak. We follow 
Ref. [53], which studied the rates in the low- and high-g 2 regions and found that for the SM 


value of eg the constraint on C 10 in Models V, VII, and VIII is 04 
(at lcr). Here it is assumed that c^/c^q M ^ is real. This implies 


< 1 + c w 


! c 


(SM) 

10 


< 1.17 


M^>y/v { zl |A 2 A 3 | x21TeV. 


(3.29) 


In the up-quark sector, LHCb established a bound Br (D + —> 1 r + /x + /x”) < 7.3 x 10~ 8 at 
90% CL [54], removing regions of m^+^- with resonance contributions. We find that this 
bound is weaker on the models considered in this paper than the one from Br(D° —>■ /i + /r _ ), 
and the latter is also theoretically cleaner. 

For FCNC top decays, CMS has set the best bound so far, Br (t —> qZ ) < 5 x 10~ 4 [55] 
at 95% CL (where q — c,u, corresponding to i — 2,1 below, respectively). Comparing to 
the dominant t —)■ bW rate, in models VI-X, 


r(t -> qZ) 2|4(! A 3 A, 
T(t bW) 


9 % 


M 4 


(3.30) 


For Models VII and VIII we used again the basis in Eqs. ( |2.32[ ) and ( J2.33| ) . We find 

M(°)/ yj'rfzl |A 3 Aj| > 0.96 TeV. The HL-LHC is expected to reach sensitivity at the 10 -5 
level [56], which will improve this bound to about 2.3 TeV. However, the direct (and A- 
independent) searches are comparably sensitive, and are expected to remain to be so. 


4■ Neutral meson mixing 

Since the new fermions interact with the Higgs field and the ordinary light quarks via 
Yukawa couplings, meson mixing, such as K°-K°, D°-D°, or B°-B 0 , is induced. The 
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FIG. 3. Diagrams contributing to K and D mixing in Model V. 


effective Hamiltonian for these processes contains dimension-six four-quark operators with 
coefficients of mass dimension — 2 . At tree level, through Z exchange, the coefficients are 
of the form ~ (AjA*) 2 u 2 /M 4 . However, coefficients of order (A,:A*) 2 /( 47 tM ) 2 are generated 
at one loop that are not CKM and/or quark-mass suppressed. For large M, these one-loop 
matching contributions are more important than tree-level Z exchange. Furthermore, they 
are independent of the Higgs vacuum expectation value, v, and arise from short distances 
~ 1/M. They can be calculated in the symmetric phase and come from box diagrams with 
virtual scalars and the heavy vector-like fermions in the loop; see Fig. [3j The resulting 
effective Lagrangians are, 


r( a ) 

^ meson 


/’(IX) 


/(X) 


/(XI) 


(a) 

-'/mix 


(X ) 2 

128vr 2 M 2 


£ HVki 7/i^jfl O {<Vrni + {47A) {dir4) 

L klmn 


+ h.c ., 


647r 2 M 2 


(A,A *) 2 

128t r 2 M 2 


-«vX) KtX) + ■ 

(urIiiUr) (XX) + h 


jxfx^*) 2 

QAtt 2 M 2 

.c., 


(X *) 2 

128t t 2 M 2 


{d/ffiAi} (XX) + h.c., 


{d l L lA) (XX) + h.c., 


(3.31) 

(3.32) 

(3.33) 

(3.34) 


with r/AA = 77 X = 1 and rj^P = = 5/4. In the Lagrangians for Models V - VIII the 

down-type quarks are the SM mass eigenstates. For the remaining models the quark fields 
occurring both in their Lagrangians and in the effective Lagrangians for meson mixing can 
be taken to be SM mass eigenstates. 

The operators above are renormalized with the subtraction point at the scale M. They 
can be related through a QCD correction factor £ to scale invariant operators whose matrix 
elements can be evaluated using lattice QCD. So we write 


{P\(q(L,R)Yq(L,R)) 2 \P) = B P f 2 P m 2 p, 


(3.35) 


where Bp does not depend on subtraction point. For K mixing, fx = 156 MeV and 
B k ~ 0.76. In that case the QCD correction factor, £, is given in the leading logarith¬ 
mic approximation wa by 


X 


r« s (M)i 

6/21 

~a s (m t y 

6/23 

'a s (m b y 

_a s (m t ) 


_a s {m b )_ 


i 

s 

Co 

O 


1 6/25 


M/«>)] 6/25 - MM)// 21 . (3.36) 


The scale does not affect any physical results. For lOTeV < M < 100 TeV, £ ~ 0.5, and 
this value is approximately the same for mixing in the K, D, and B ds systems. We use the 






















17 


lattice QCD averages from Ref. [50], and B D = R^ s (3 GeV)/[a s (3 GeV )] 6 / 25 ~ 1.0 from a 
recent calculation of R^ s (3GeV) [58] . 

In the case of kaon mixing, we demand the new physics contribution to the real part of 
M 12 to be less than 40% of the experimental value of A mx [59] • We take the la uncertainty 
on the determination of ex in the SM to be 25% (reading off the bound hx < 0.35 from 
Ref. [ 6 U]). For all the (hadronic) models, except Model X, this results in the constraint 


M (a) 


> max 


42 


>Jv mil |Re(A 2 A* 1 ) 2 | , 6.7 x 10 2 \Jr ££ |lm(A 2 A ];) 5 


TeV. 


(3.37) 


(Note that the 90% CL constraint from ex is only slightly weaker, replacing 6.7 x 10 2 TeV by 
6.3 x 10 2 TeV.) In the future, the ex constraint is e xpect ed to improve by about \[ 7 1A [60]. 
replacing 6.7 x 10 2 TeV by 1.0 x 10 3 TeV in Eq. (3.37). Improvement in the real part 


is contingent upon lattice QCD calculations of the long-distance contributions to A m K ', 
reaching x% precision would replace 42 TeV in Eq. (3.37) by 260 TeV/ y/x. 

In B ds meson mixing, the new physics contribution is conventionally parametrized as 
M ±2 = (IW^Qsm x (l + h dtS e 2lcrd ’ s ). Until recently, the bounds on real (2a = 0, mod7r, that 
is MFV-like) and imaginary (2a = n/2, mod7r) new physics contributions have been quite 
different [HDE2]. This is no longer the case [60], and since we are most interested in physics 
reach, we simply quote the limits on the absolute values of the new physics contribution. 
In B d ,s mixing, h d < 0.3 and h s < 0.2 [60] yield in all models except Model X, using 
f Bd = 188MeV, f Ba = 226MeV, B Bd = 1.27, B Bs = 1.33 [50], 


M< a) 


> 


it|A 3 Ai| 2 x25TeV, 


AfW 


> 


VmL IA 3 A 2 I 2 x 6.4 TeV. 


(3.38) 


(For the imaginary part in B d mixing, 25 TeV should be replaced by 31 TeV, which is a 
smaller dependence on the phase of new physics than those ignored for A F = 1 FCNC 
transitions earlier.) In the above constraints = 2, = 1, and for Model IX it is the 

y«P s ^hat occur j n ^} ie constraints. I 11 the next decade these limits will improve to h d < 0.05 


and h s < 0.04 [60], which will replace 25 TeV by 61 TeV and 6.4 TeV by 14 TeV in Eq. (3.38). 

The mixing of D mesons is probably dominated by long-distance physics. Thus, we 
can only require that the new physics contribution does not exceed the measurement, i.e., 
Am D /Y D < 0.006 [3TJ. (The significance of A m B % 0 is less than 2a, so we use the upper 
bound of the la region. We do not distinguish between imaginary and real contributions to 
M 12 relative to the SM; the bounds on the mass scale may differ by a factor ~ 2, depending 
on the value of A m B .) In contrast to above, the A^ couplings occur in the constraint 
for Model IX. Furthermore, in this case we choose the up-type quarks to be the SM mass 
eigenstate fields, so the factors of the CKM matrix move to the terms with the down-type 
quarks. Then D mixing implies for all models except Model XI, 


M (a) 


> 


^mii |A 2 A;l | 2 x 48 TeV, 


(3.39) 


CX.) 

where now r/ ini (. = 1. Note that for Models V-VIII the constraint from D mixing is slightly 
stronger than from the real part of K mixing. However, the constraints are actually a little 
different since the A’s in Eqs. ( |3.38 ) and in (3.39) are not the same. They are linearly related 
through the CKM matrix. The future evolution of this bound is uncertain. While Belle II 
expects to measure Am^/T B with an uncertainty of 0.001 [28] . the central value will matter 
for the bounds on new physics, and therefore we do not assume that this bound will improve. 
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5. Unitarity of the CKM Matrix 


The CKM matrix is not unitary in the extensions of the SM we are considering. The 
strongest constraints on violations of unitarity come form the first row and first column of 
the CKM matrix m 

Sn = 0.9999 ± 0.0006 and T u = 1.000 ± 0.004. (3.40) 


In the models with vector-like quarks both of these constraints involve the same combination 
of couplings and we find 

M (a) > |Ai| \frff x 7.8 TeV , (3.41) 

where in = 277un H) = = 1 and rjun' = — 0- For the second 

row we have the constraint 

S 22 = 1.002 ±0.027, (3.42) 


which yields 

MW > |A 2 |y^X 1.1 TeV. 


(3.43) 


IV. CONCLUSIONS 

There are 11 renormalizable models that add to the SM vector-like fermions in a single 
(complex) representation of the gauge group that can Yukawa couple to the SM fermions 
through the Higgs field. These BSM fermions can have a mass M that is much greater 
than the weak scale, since they have a mass term even in the absence of weak symmetry 
breaking. However, unlike BSM scalars, such fermions are technically natural. These models 
are a class of very simple extensions of the SM that do not worsen the SM hierarchy puzzle. 

The masses of these vector-like fermions can take any value, up to the ultraviolet cutoff 
(~ Mpi), and so there is no particular reason that they should be in a region that can be 
probed experimentally. However, there are many such models and it is not unreasonable 
that one of them has vector-like fermions with masses in the experimentally testable range. 

We considered the experimental constraints from flavor physics on the mass of these 
vector-like fermions to get a feel for the mass reach that the present experiments have for 
this class of models. We are primarily interested in very heavy vector-like fermions, say, 
with masses greater than 10 TeV. Hence it is flavor-changing neutral-current processes that 
provide the most important constraints. However, we also discussed violations of lepton 
universality and CKM matrix unitarity. 

An important feature of these models is that for large M the BSM contribution to meson 
mixing is either suppressed by a loop factor (~ 1/167T 2 ) or a factor of m 2 z /M 2 compared 
to processes that change flavor by one unit and are dominated by tree level through flavor¬ 
changing Z exchange. This implies that, except for the case of the kaon CP violation 
parameter ex, the constraints on M from meson mixing are not overwhelmingly strong. We 
computed the order 1/(47 tM) 2 one-loop contribution to the coefficients of the four-quark 
operators responsible for meson mixing that is not suppressed by SM quark masses or weak 
mixing angles. 
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Model Quantum 

Present bounds on Mj TeV and XiXj for each ij pair 

numbers 

ij = 12 

ij = 13 

ij = 23 

I (1,1,-1) 

310“ 

7.& 

7.4 C 

II (1,3, —1) 

220“ 

4.9 fc 

5.2 C 

III (1,2,-1/2) 

310“ 

7.0 fc 

7.4 C 

IV (1,2,-3/2) 

310“ 

7.0 6 

7.4 C 


AF = 1 AF = 2 

AF = 1 AF = 2 

AF = 1 AF = 2 

V (3,1,-1/3) 

m d [100] e (42, 670}7 

30^ 25 /l 

21* 6.43 

VI (3,1,2/3) 

3.9 fc (42, 670}7 

2 h h 

6.43 

VII (3, 3,-1/3) 

47 d [71] e {47, 750}7 

2R 28 /l 

15* 7.23 

VIII (3, 3, 2/3) 

66 d [100] e {47, 750}7 

303 28 /l 

21* 7.23 

IX A(“) T / 

3.9 fc 67 z 


9.13 

IX AC) }^’ i/b) 

66 d [100] e {59, 950}7 

303 35 h 

18 m 9.13 

X (3,2,7/6) 

3.9 fc 48* 

— — 

— — 

XI (3, 2-5/6) 

66 d [100] e {42, 670}7 

303 25 h 

18 m 6.43 


TABLE I. Bounds from flavor-changing neutral currents on M [TeV]/-y/|AjAjl in the leptonic mod¬ 
els, and from the A F = 1 constraints on the hadronic models. The A F = 2 bounds show 
M /y 7 )A ? ;Aj| 2 ; except for K meson mixing we show { M/ y // |Re(AjAj) 2 |, M /y/|InL(AjAj) 2 | }. The 
strongest bounds arise from: a) fi to e conversion; b) r — > e7r; c) r —>• pp; d ) K — > 7rz/u; e) 
AT —> p + (this bound involves |Re(AiA 2 )|); f) K mixing; g) B —>• irp + p~-, /i) B d mixing; ?') 
B -» X s £ + £ _ ; j) B s mixing; &) F —>• p + l ) F mixing; m) B s —> p + g~. 


We are interested in a rough assessment of the experimental reach. The strongest current 
bounds on the vector-like fermion masses and couplings in each of the 11 models studied in 
this paper are summarized in Table [Tj We display bounds which are above or will get near the 
10 TeV level in the near future. Muon to electron conversion in a nucleus, fi —> 3e, the kaon 
CP -violating parameter 6k, and Kl -A pC are sensitive to vector-like fermion masses ^ 
100 TeV (for Yukawa couplings with magnitude unity). However, it is important to remember 
that the couplings Aj may have a flavor structure that suppresses these contributions relative 
to those involving the third generation. 

We summarize the expected future sensitivity in Table [TTJ where we only display sensitiv¬ 
ities near or above 10 TeV. The Mu2e constraint will be improved dramatically in the next 
generation of experiments [221123]. The measured K + —y 7 Cuu branching ratio corresponds 
to a mass reach around 70 TeV, which will increase substantially as the next generation 
experiments reach an uncertainty at about 10% of the SM rate, especially since the current 
central value is above the SM prediction (which has very small theoretical uncertainty). 
The improvement in the A m,x bound is entirely dependent on lattice QCD calculations, as 
discussed after Eq. (3.37). The B ds and D mixing sensitivities will be improved by Belle 11 
and LHCb. These experiments, and CMS and ATLAS, will also probe FCNC B d s , D } and 
r decays much better than current bounds. The future sensitivities in Table [TT] correspond 
to estimated 50/ab Belle If [28] and 50/fb LHCb ^5] sensitivities and CMS/ATLAS reach 
in rare decays on the same time scale. Compared to Table [IJ a greater number of the best 
bounds will come from purely leptonic rather than semileptonic decays. The sensitivities in 
Table [TT] may be realized in ~ 10 years. 
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Model 

Quantum 


Future bounds on M/TeV and A* A j 

for each ij pair 


numbers 


ij = 12 

ij = 13 


ij = 23 


I 

(1,1,-1) 


2000“ 


19 b 



21 c 


II 

(1,3, -1) 


1400“ 


13 b 



15 c 


III 

(1,2,-1/2) 


2000“ 


19 b 



21 c 


IV 

(1,2,-3/2) 


2000“ 


19 b 



21 c 




A F = 

1 AF = 2 

AF = 

1 AF = 

2 

A F = 

1 AF = 2 

V 

(3,1,-1/3) 

280 d 

{100, 1000} e 

607 

613 


39 h 


14* 

VI 

(3,1,2/3) 

8..P 

{100, 1000} e 

— 

613 


— 


14* 

VII 

(3, 3,-1/3) 

200 d 

{110, 1100} e 

427 

683 


28 h 


16* 

VIII 

(3, 3, 2/3) 

280 d 

{110, 1100} e 

6(]7 

683 


39 h 


16* 

IX \V) 

| (3, 2,1/6) 

8.:v 

67 k 

— 

863 


— 


20* 

IX A( rf ) 

280 d 

{140, 1400} e 

607 

863 


39 h 


20* 

X 

(3, 2, 7/6) 

8.3 j 

48 fc 

— 

— 


— 


— 

XI 

(3, 2-5/6) 

280 d 

{100, 1000} e 

607 

613 


39 h 


14* 


TABLE II. Expected future bounds from flavor-changing neutral currents on M [TeV]/-y/| \i\j\ in 
the leptonic models, and from the A F = 1 constraints on the hadronic models. The A F = 2 bounds 
show M /-^/lAjAjl 2 , except for K meson mixing we show { M/ y^|Re(A,:A*) 2 |, M /y/|Im(AjAj) 2 | }. 
The bounds are from: a) n to e conversion; b) r —> e7r; c) r —> fip: d ) K —>• nisi?-, e ) K mixing; /) 
B d g + g~; g) B d mixing; h) B s g + g~] i ) B s mixing; j ) D ->• fe) D mixing. 
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Appendix A: Diagonalizing matrix 


Here we summarize the 4x4 diagonalizing matrices in the 11 models. Terms of order 
( v/M) 2 and higher are not explicitly displayed. 


Model I : VI 
Model II : Vl 


1 A i.v/M 
—A* v/M l 3 x 3 


' H — 1-4x4 , 


1 -\ iV /V2M \ 

A *v/V2M l 3 x 3 )’ ^- l4x4 ’ 


(Al) 

(A2) 



Model III : VI 
Model IV : V[ 
Model V : V[ 
Model VI : V}'; 

Model VII : V}'; 

vt 

Model VIII : V“ 

vi 

Model IX : V£ 

vi 

Model X : V£ 
Model XI : v[ 


14x4 j Vr — 


14x4 , Vr ~ 


1 -A *v/M 

A iv/M l 3x3 

1 A *v/M 

—A i,v/M l 3x3 

A ^/M 


-A *v/M l 3x3 

1 —A iv/M 
A *v/M l 3x3 

1 A iv/M 

—X*v/M l 3x3 


, V2 = 1 4x4 , 

j Vr = 14x4 ) 
, Vr — 1 4X 4 ; 


( 1 -\v/V2M\ 

V A *v/V2M l 3x3 ) ’ ' 


= 1 


R — J-4X4 ? 


, V* = I 4x4 , 


1 -A iv/y/2M \ 

\*v/V2M l 3x3 J ’ 

1 —A iv/M 
A *v/M l 3x3 

1 A \ u) *v/M 
-A [ U) v/M l 3x3 

1 -A [ d) *v/M 


Vr — 14x4 5 


Ux4 5 v R 


v% = 


14x4 : 


V d - 

V R ~ 


L 4x4 i Vr 


14x4 : 


v u - 

v R — 
V d — 

V R ~ 


\i ] v/M l 3x3 

1 -A *v/M 

X iv/M l 3x3 

1 —X*v/M 
Xiv/M l 3x3 
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(A3) 

(A4) 

(AS) 

(A 6 ) 

(A7) 

(AS) 

(A9) 

(A 10 ) 

(AH) 

(A 12 ) 

(A13) 

(A14) 


Finally, there are tree-level mass splittings among the heavy fermions, except for the 
SU(2)r singlet models, due to electroweak symmetry breaking. The results are 


Model II : M 0 = M + AM , M_ x — M + AM/2 , M _ 2 = M , (A15) 

Model III : M 0 = M , M_ x = M + AM , (A16) 

Model IV : = M + AM , M _ 2 = M, (A17) 

Model VII : M 2/3 = M + AM , M _ 1/3 = M + AM/2 , M _ 4/3 = M , (A18) 

Model VIII : M 5/3 = M , M 2/3 = M + AM/2 , M _ 1/3 = M + AM , (A19) 

Model IX : M 2/3 = M_ 1/3 = M + AM , (A20) 

Model X : M 5/3 = M, M 2/3 = M + AM, (A21) 

Model XI : M _ 1/3 = M + AM, M _ 4/3 = M. (A22) 


In each model AM = (v 2 /2M) |A/ 2 
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Appendix B: Numerical inputs and Wilson coefficients 


We collect in this Appendix expressions used to derive bounds on vector-like fermions 
from flavor-changing neutral-current processes, to facilitate easier comparison with and re¬ 
production of our numerical results. Our goal in this paper is to study the sensitivities of 
many processes, so leading or next-to leading order results suffice. (For most processes the 
state of the art is one or two orders higher.) In many cases we ignore the SM uncertainties, 
when we know that they are subdominant effects. 

We adopt for the numerical values of the coupling constants a:(m&) = 1/133, sin 2 9w = 
0.23, gz = 0.73. For the top quark mass we use m t (m t ) = 165 GeV, obtained from the 
one-loop relation from the m t = 173 GeV (presumed) pole mass, extracted from fits to tt 
production at the Tevatron and the LHC. 

The SM Wilson coefficients are as follows. For K + —y 7r + i , i 9 i [361 HU] 


(SM) _ (SM) 
c +v ~ G — 


V2G 


F OL 


v t :v ld x(x ,)+ 


7 r sin 9 W l 

where x t = fn 2 /rri 2 Vl A = 0.225, P C (X) — 0.4 [38| HO] F] and 


x(x t ) = x 4 


x t + 2 (3x t - 6 ) In xt 
_x t - 1 (x t - l ) 2 


(Bl) 


(B2) 


Here we have neglected the electroweak corrections. 
For B q —* ji + n~ (where q = s,d), 


jsu) _ r F 


G F OL 


C 10 — 


Y{x t ) 

a/2 7T tb tq sin 2 6 W 


V*V tn 


Y(x t ) = 


Xt ( x t — 4 3ay In x t 


+ 


x t - 1 (x t - 1 y 


(B3) 


(In the usual notation, Cio = —Y(x t )/ sin 2 9w — —4.2.) 

The short-distance contribution to the Kl —>• fi + / jl~ rate can be written as 


r (K l —>■ g + g, )sd — 


G i- 4 


F m w ■ 4 Q tl 


47T 5 


sin 9 w fj < m K m 


Arri'i 


Re 


m 


K 


v*v td Y(x t ) + v; s v cd \ 4 p c 


2 1 


where P c ~ 0.11. 

The B —> rate is given by 021 


(B4) 


dr(B -> A>'/» 

dg 2 


2 

Lt r 


? a 7TV 


B 


15367T 5 
4m 2 


\VM,\ 2 [A K (, 2 )] 3 /q Sl [|Cf(9 2 )| 2 + |C W + c; 0 | 2 ] 

4 nib 


(m B + m K ) 2 


fr l^rl + 


m B + m K 


M+Re[Cf( 9 2 )C 7 *] (B5) 


where A k(q 2 ) is a phase space factor, f + and are g 2 -dependent form factors. The B —> 
7 r fjpn~ rate is obtained with obvious replacements. In the heavy quark limit ./'+ //t = 
1 + 0 (AQCD/iTib ) [64], and model calculations are consistent with a mild q 2 dependence of 


8 In our definition of c/ SM '* all lepton flavors (labeled by i) are included by using parameter P C (X), which 
is defined in ES1I33- Then c/ SM ^ corresponds to the decay amplitude. 
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this ratio. Hence, at the desired level of precision, the form factors can be pulled out of the 
expression in Eq. (B5), obtaining simple approximations for the effect of new physics 
via Cio and C' w . (This is impossible for B —* K*/d + fx~, as CV enters with a 1/q 2 dependence 
in that case.) In the numerical analysis, we use C 7 = —0.33, Cio = —4.2, and a mean value 
|Cf(g 2 )| = 4.4. These estimates can be refined as the measurements improve. 


[1] G. C. Branco and L. Lavoura, Nucl. Phys. B 278, 738 (1986). 

[2] Y. Nir and D. J. Silverman, Phys. Rev. D 42, 1477 (1990). 

[3] F. del Aguila, L. Ametller, G. L. Kane and J. Vidal, Nucl. Phys. B 334, 1 (1990). 

[4] G. Barenboim, F. J. Botella and O. Vives, Nucl. Phys. B 613, 285 (2001) [hep-ph/0105306]. 

[5] G. Cacciapaglia, A. Deandrea, D. Harada and Y. Okada, JHEP 1011, 159 (2010) 
[arXiv: 1007.2933 [hep-ph]]. 

[6] W. Altmannshofer, M. Bauer and M. Carena, JHEP 1401, 060 (2014) [arXiv: 1308.1987 [hep- 
ph]]. 

[7] T. Ma, B. Zhang and G. Cacciapaglia, Phys. Rev. D 89, no. 1, 015020 (2014) [arXiv: 1309.7396 
[hep-ph]]. 

[8] S. A. R. Ellis, R. M. Godbole, S. Gopalakrishna and J. D. Wells, JHEP 1409, 130 (2014) 
[arXiv: 1404.4398 [hep-ph]]. 

[9] F. del Aguila, M. Perez-Victoria and J. Santiago, JHEP 0009, 011 (2000) [hep-ph/0007316]. 

[10] A. K. Alok, S. Banerjee, D. Kumar, S. U. Sankar and D. London, arXiv:1504.00517 [hep-ph]. 

[11] G. Aad et al. [ATLAS Collaboration], arXiv: 1503.05425 [hep-ex], 

[12] A. Succurro [ATLAS Collaboration], EPJ Web Conf. 60, 20037 (2013). 

[13] V. Khachatryan et al. [CMS Collaboration], arXiv: 1503.01952 [hep-ex]. 

[14] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 112, no. 17, 171801 (2014) 
[arXiv: 1312.2391 [hep-ex]]. 

[15] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 729, 149 (2014) [arXiv:1311.7667 
[hep-ex]]. 

[16] G. Barenboim and F. J. Botella, Phys. Lett. B 433, 385 (1998) [hep-ph/9708209]. 

[17] K. Ishiwata and M. B. Wise, Phys. Rev. D 88, no. 5, 055009 (2013) [arXiv: 1307.1112 [hep-ph]]. 

[18] F. J. Botella, G. C. Branco and M. Nebot, JHEP 1212, 040 (2012) [arXiv: 1207.4440 [hep-ph]]. 

[19] S. Fajfer, A. Greljo, J. F. Kamenik and I. Mustac, JHEP 1307, 155 (2013) [arXiv:1304.4219 
[hep-ph]]. 

[20] G. Cacciapaglia, A. Deandrea, N. Gaur, D. Harada, Y. Okada and L. Panizzi, arXiv: 1502.00370 
[hep-ph]. 

[21] W. H. Bertl et al. [SINDRUM II Collaboration], Eur. Phys. J. C 47, 337 (2006). 

[22] R. J. Abrams et al. [Mu2e Collaboration], arXiv:1211.7019 [physics.ins-det]. 

[23] A. Kurup [COMET Collaboration], Nucl. Phys. Proc. Suppl. 218, 38 (2011). 

[24] A. de Gouvea and P. Vogel, Prog. Part. Nucl. Phys. 71, 75 (2013) [arXiv: 1303.4097 [hep-ph]]. 

[25] U. Bellgardt et al. [SINDRUM Collaboration], Nucl. Phys. B 299, 1 (1988). 

[26] K. Hayasaka et al. [Belle Collaboration], Phys. Lett. B 687, 139 (2010) [arXivTOO 1.3221 
[hep-ex]]. 

[27] A. Bravar [Mu3e Collaboration], Nucl. Part. Phys. Proc. 260, 155 (2015). 

[28] B. Golob, K. Trabelsi, P. Urquijo, BELLE2-NOTE-0021, https://belle2.cc.kek.jp/ 
~twiki/pub/B2TiP/WebHome/belle2-note-0021.pdf. 





24 


[29] Y. Miyazaki et al. [Belle Collaboration], Phys. Lett. B 672, 317 (2009) [arXiv:0810.3519 [hep- 
ex]]. 

[30] K. Hayasaka, arXiv: 1010.3746 [hep-ex]. 

[31] Y. Amhis et al. [Heavy Flavor Averaging Group Collaboration], arXiv:1412.7515 [hep-ex]; and 
updates at http://www.slac.stanford.edu/xorg/hfag/, 

[32] A. Aguilar-Arevalo et al. arXiv: 1506.05845 [hep-ex], 

[33] S. Ito et al, J. Phys. Conf. Ser. 631, no. 1, 012044 (2015). 

[34] C. Lazzeroni et al. [NA62 Collaboration], Phys. Lett. B 698, 105 (2011) [arXiv: 1101.4805 
[hep-ex]]. 

[35] A. J. Buras, F. De Fazio and J. Girrbach, JHEP 1302, 116 (2013) [arXiv:1211.1896 [hep-pli]]. 

[36] G. Buchalla and A. J. Buras, Nucl. Phys. B 412, 106 (1994) [hep-ph/9308272], 

[37] G. Buchalla and A. J. Buras, Nucl. Phys. B 548, 309 (1999) [hep-ph/9901288]. 

[38] A. J. Buras, F. Schwab and S. Uhlig, Rev. Mod. Phys. 80, 965 (2008) [hep-ph/0405132v3]. 

[39] J. Brod and M. Gorbahn, Phys. Rev. D 78, 034006 (2008) [arXiv:0805.4119 [hep-ph]]. 

[40] A. J. Buras, D. Buttazzo, J. Girrbach-Noe and R. Knegjens, arXiv: 1503.02693 [hep-ph], 

[41] A. V. Artamonov et al. [E949 Collaboration], “New measurement of the K + — > ir + vv branch¬ 
ing ratio,” Phys. Rev. Lett. 101, 191802 (2008) [arXiv:0808.2459]. 

[42] G. Buchalla, G. Hiller and G. Isidori, Phys. Rev. D 63, 014015 (2000) [hep-ph/0006136]. 

[43] W. Altmannshofer, A. J. Buras, D. M. Straub and M. Wick, JHEP 0904, 022 (2009) 
[arXiv:0902.0160 [hep-ph]]. 

[44] V. Khachatryan et al. [CMS and LHCb Collaborations], arXiv:1411.4413 [hep-ex], 

[45] LHCb Collaboration, LHCb-PUB-2014-040, https://cdsweb.cern.ch/record/1748643, 

[46] M-O. Bettler, presentation at “Physics at the High-Luminosity LHC”, https://indico. 
cern.ch/event/360104/, 

[47] K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 (2014). 

[48] G. Isidori and R. Unterdorfer, JHEP 0401, 009 (2004) [hep-ph/0311084]. 

[49] R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 725, 15 (2013) [arXiv: 1305.5059 [hep-ex]]. 

[50] S. Aoki et al. [Flavor Lattice Averaging Group], Eur. Phys. J. C 74, 2890 (2014) 
[arXiv:1310.8555 [hep-lat]]. 

[51] Tim Gershon, private communications. 

[52] R. Aaij et al. [LHCb Collaboration], JHEP 1212, 125 (2012) [arXiv:1210.2645 [hep-ex]]. 

[53] T. Huber, T. Hurth and E. Lunghi, arXiv: 1503.04849 [hep-ph], 

[54] R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 724, 203 (2013) [arXiv: 1304.6365 [hep-ex]]. 

[55] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 112, no. 17, 171802 (2014) 
[arXiv: 1312.4194 [hep-ex]]. 

[56] ATLAS Collaboration, arXiv: 1307.7292 [hep-ex], 

[57] F. J. Gilman and M. B. Wise, Phys. Lett. B 93, 129 (1980). 

[58] N. Carrasco et al. [ETM Collaboration], arXiv:1505.06639 [hep-lat], 

[59] Z. Bai, N. H. Christ, T. Izubuchi, C. T. Sachrajda, A. Soni and J. Yu, Phys. Rev. Lett. 113, 
112003 (2014) [arXiv: 1406.0916 [hep-lat]]. 

[60] J. Charles, S. Descotes-Genon, Z. Ligeti, S. Monteil, M. Papucci and K. Trabelsi, Phys. Rev. 
D 89, no. 3, 033016 (2014) [arXiv: 1309.2293 [hep-ph]]. 

[61] K. Blum, Y. Grossman, Y. Nir and G. Perez, Phys. Rev. Lett. 102, 211802 (2009) 
[arXiv:0903.2118 [hep-ph]]. 

[62] G. Isidori, Y. Nir and G. Perez, Ann. Rev. Nucl. Part. Sci. 60, 355 (2010) [arXiv: 1002.0900 
[hep-ph]]. 



25 


[63] A. J. Buras, F. De Fazio, J. Girrbach, R. Knegjens and M. Nagai, JHEP 1306, 111 (2013) 
[arXiv: 1303.3723 [hep-ph]]. 

[64] N. Isgur and M. B. Wise, Phys. Rev. D 42, 2388 (1990). 



